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Abstract. We describe a work in progress on the reconstruction of radial profiles 
for the ion effective charge Z^ff on the TEXTOR tokamak, using statistical 
data analysis techniques. We introduce our diagnostic for the measurement 
of bremsstrahlung emissivity signals. Z^s profiles can be determined by Abel- 
inversion of line-integrated measurements of the bremsstrahlung emissivity (efj) 
from the plasma and the plasma electron density (rie) and temperature (To). 
Ifowever, at the plasma edge only estimated values are routinely used for ric and 
Te, which are moreover determined at different toroidal locations. These various 
uncertainties hinder the interpretation of a Z^g profile outside the central plasma. 
In order to circumvent this problem, we propose several scenarios meant to allow 
the extraction by (Bayesian) Blind Source Separation techniques of either (line- 
integrated) Z(,ff waveshapes or absolutely calibrated signals from (line-integrated) 
emissivity signals, using also density and temperature signals, as required. 



1. The role of impurities and Zcft in a tokamak plasma 

The behavior of plasma impurities in a tokamak plasma is a critical issue On 
the one hand, these impurities increase the fuel dilution and, at reactor parameters, 
they can be responsible for a considerable power loss from the plasma core by 
bremsstrahlung. On the other hand, line radiation by \ow-Z impurities at the plasma 
edge can be advantageous for the reduction of the plasma-wall interaction. 
To study the behavior of impurities in the plasma, one needs a local measure for 
impurity concentration. A parameter of great interest is the so-called ion effective 
charge of the plasma, Z^s, which is averaged over all species. Once the plasma electron 
density (rie) and temperature (Te) are known, it is proportional to the bremsstrahlung 
emissivity (eg) from the plasma 0, which originates mainly from the acceleration of 
electrons in the field of ions: 



and where Ui and Zi are the density and charge of ion species i. 
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2. A diagnostic for the determination of Zcs on TEXTOR 



At the tokamak TEXTOR (FZJ Jiilich, Germany), we run a diagnostic for the 
determination of ZeS from measurements of bremsstrahlung emissivity in the visible 
(see also 0). 21 fiber optic cables are directed along different lines of sight in a 




Figure 1. The viewing geometry for the bremsstrahlung Z^ff diagnostic at 
TEXTOR. The dashed lines represent two of the viewing chords used in the ICA 
analysis described below. 



poloidal cross-section of the plasma (figure^]). The hght is led through an interference 
filter (peak transmission at 523 nm) and focused onto the CCD of a cooled 12-bit 
camera (DTA Inc., Italy). The transmitted wavelength region is known to consist for 
99% of bremsstrahlung continuum. The measured emissivity signals are stored into 
the internet-accessible TEXTOR central database. 

In addition, four other channels will be equipped with a photomultiplier, in order to 
study some phenomena taking place on smaller time scales as well. We expect to have 
this second system operational by the beginning of the TEXTOR campaign starting 
in the last quarter of 2004. 

A central line-averaged ZeS is routinely computed and is available online. The average 
Zcs is defined as (see 
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where a is the plasma radius and the so-called (averaged) Gaunt factor. In 
addition, a radial profile for Z^ff (figure ^ can be reconstructed from the full set 
of line- integrated emissivities by an Abel- inversion [2] . This also requires profiles for 
electron density and temperature. 



3. Difficulties in the reconstruction of Z^ff profiles 

Several factors render the radial matching of the different profiles problematic, 
hindering the interpretation of a Zcff profile outside the central plasma. Indeed, at 
the plasma edge only estimated values are routinely used for ric and To, which are 
moreover determined at different toroidal locations. As a consequence, one frequently 
observes an unphysical rise of Z^ff toward the plasma edge (figure [S)). In fact, so far 
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Figure 2. A typical Z^ff profile obtained by Abel-inversion of line-integrated 
quantities. The error bars are obtained from the estimated errors on the 
emissivity, density and temperature profiles, as described in (S). 
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Figure 3. An extreme case of the problems that might arise in the interpretation 
of a reconstructed Z^g profile. The profile diverges at the inboard side of the 
tokamak. 



none of the available methods for the determination of Z^s has provided a full profile, 
which is at present a real challenge. 

Conversely, if it were e.g. possible to obtain a set of line-integrated values for Zcs 
directly from the line-integrated measurements of es, no and Te, then these problems 
would be avoided, as a ZcS profile could directly be obtained from inversion of the 
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Figure 4. A line-integrated emissivity signal (solid curve) and its fit using a 
linear regression of a line-integrated eg on line-integrals (same line of sight) of 
and Zeff . Despite of the several simplifications, the regression model seems to 
hold. 



line-integrals. Alternatively, one may work directly with already inverted data. As 
explained below, this even opens the possibility to obtain estimations for Zcs without 
using density and temperature information. 

In a first attempt toward these goals, we want to get an idea of the relative signal of 
ZeB using bremsstrahlung emissivity measurements and if necessary also density and 
temperature signals. 

4. Extraction of a Zcff signal from line-integrated bremsstrahlung 
measurements; test using linear multi-channel ICA 

From equation , the signal for Z^s can be seen as a nonlinear mixture of the signals 
for eg, no and To- The dependence on Tg is, however, weak and the To-dependence can 
in a first approximation be left out of the analysis. We may linearize the expression 
of eg as a function of and Z^s up to first order, i.e.: 

Eff sa Ci?ic + C-iZcff, Ci, C2 regression coefficients, (3) 

which constitutes a linearized mixing model. The traditionally obtained signals for the 
values of the plasma parameters can be used to perform a linear regression analysis, 
and a similar analysis can be performed for the corresponding line-integrated signals, 
upon radial integration of the model The results are displayed in figure^ 
As a preliminary test (see also U and [5]) for the process of extraction of a Z^s signal 
from bremsstrahlung emissivity signals, four emissivity line-integrals along different 
lines of sight were fed to the FastICA program |B] for linear Independent Component 
Analysis (ICA) [T]. 

I In fact, one finds that one may even use directly ne (and Tc) as input signals for the regression. 



Reconstruction of Z^s profiles at TEXTOR through Bayesian source separation. 

I 



4 



^ - 

v. 



V-. 



■ ell ^ 



in 1 

Time (<\ 



Figure 5. The time trace of the hne-integrated emissivity from a central 
channel (ch 4) and a more peripheral channel (eh 5) (shown in figure 0. The 
signal (rescaled) of the central line-integrated density is also shown. Two striking 
features in the signal are marked by dashed lines. 



ICA is a technique that can be derived in a Bayesian context from first principles. In 
(Hnear) ICA, it is assumed that a number n of signals Xi{t) is measured, which are 
each a linear mixture of (ideally) the same number of unknown source signals Si{t) 
that are assumed to be statistically independent: 
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where A is an unknown mixing matrix. ICA now consists of estimating both the 
source signals and the mixing matrix. The idea is in principle to find a set of linear 
combinations of the xt that are maximally nongaussian. Due to the central limit 
theorem these will equal the independent components Si. Thus, the ICA procedure 
requires the maximization of some measure of nongaussianity (e.g. negentropy), which 
is accomplished by some gradient descent approach. 

In our application, it is thus implicitly assumed that the waveform for any time trace 
is the same at different radii, which is on the average true, and any local diversion 
from this average will be considered as due to noise, or is neglected while taking 
different regression coefficients for different lines of sight. The results of this analysis 
are displayed in figures [S] and (HJ and can be compared to the central signal for ZcS 
(figure El . In this analysis, no density or temperature information was used. 
It should be noted however that the signals for the plasma parameters are, to a certain 
extent, interrelated, and thus the assumption of independence of the components is 
not entirely fulfilled in this instance. Also, the obtained Z^s waveform is in fact an 
average over a certain viewing angle in the poloidal cross-section, and can thus hardly 
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Figure 6. The time trace (ICA) of one of the extracted independent components 
and tlie signal (rescaled) of the central Z^ff. Note the good correspondence, 
especially at the time of the bumps in the emissivity signal. 

be used for the reconstruction of a reliable ZcS profile. These issues are treated in the 
next section. 

5. Outlook: Bayesian source separation models 

The signals for eg, Uo, To and Z^ff are mutually dependent, so the independence 
assumption of the ICA scheme does not hold, although in practice it still yields 
acceptable results, as demonstrated above. We would however like to have more 
control over the actual source signal that is extracted from the mixture, in order to 
obtain better estimations of a Zcs waveform during discharges in different plasma 
regimes. 

Here, the Bayesian line of thought might come to help, as it allows us to incorporate 
prior information on the physics of the problem and the ZoS waveshape into the 
separation model. This information can be extracted from the central and line- 
averaged ZcS signals obtained via traditional methods. 

In practice, the aim in Bayesian source separation is to obtain the posterior distribution 
of the unknown signals S in terms of the likelihood of the measured data X and the 
prior information (encoded in p{S)) one is willing to assume: 

p{S\X)^piX\S)p{S). (5) 

We are thus looking for signals S that maximize this posterior probability, given the 
prior distribution and the measured data. We then also use known theoretical relations 
between the data and the signals S, such as 

This analysis can be conducted as a single-channel method (see :8 ) on line-integrated 
signals, and when it is done for several lines of sight, it allows us to directly obtain 
a ZcS profile via Abel-inversion. As mentioned before, instead one can also operate 
immediately on profile data. 
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6. Conclusion 

Several difficulties are associated with the calculation of Zes profiles in a tokamak 
plasma, in particular at the plasma edge. Determining ZeS signals directly from the 
brcmsstrahlung data can avoid these problems. We have conducted a preliminary test 
with good results for the construction of Zgg line-integrals via Independent Component 
Analysis. Finally, we have argued that the separation process can be made more robust 
using a Bayesian source separation technique, which can employ prior information on 
the Zgg signal. 
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